18.7 KB · ⇊ 6 parts · ⤓ Micron 18.7 KB · ⤓ HTML 171.6 KB

top

Die Itō-Formel (auch Itō-Döblin-Formel; selten auch Lemma von Itō), benannt nach dem japanischen Mathematiker Itō Kiyoshi, ist eine zentrale Aussage in der stochastischen Analysis. In seiner einfachsten Form ist es eine Integraldarstellung für stochastische Prozesse, die Funktionen eines Wiener-Prozesses sind. Es entspricht damit der Kettenregel bzw. Substitutionsregel der klassischen Differential- und Integralrechnung.

Itô publizierte 1951 einen Beweis.cite-ref-1[1]

• Version für Wiener-Prozesse

• Version für Itō-Prozesse

• Mehrdimensionale Version

• Version für Semimartingale

• Bemerkung

• Für das Stratonowitsch-Integral

• Version für Funktionen mit beschränkter quadratischer Variation

• Beispiele

• Unendlich-dimensionale Itō-Formeln

• Siehe auch

• Literatur

• Einzelnachweise

Sei ( W t ) t ≥ ≥ 0 {\displaystyle (W_{t})_{t\geq 0}} ein (Standard-)Wiener-Prozess und h : : R → → R {\displaystyle h\colon \mathbb {R} \to \mathbb {R} } eine zweimal stetig differenzierbare Funktion. Dann gilt

h ( W t ) = h ( W 0 ) + ∫ ∫ 0 t h ′ ( W s ) d W s + 1 2 ∫ ∫ 0 t h ″ ( W s ) d s . {\displaystyle h(W_{t})=h(W_{0})+\int _{0}^{t}h'(W_{s})\,{\rm {d}}W_{s}+{\frac {1}{2}}\int _{0}^{t}h''(W_{s})\,{\rm {d}}s\,.}

Dabei ist das erste Integral als Itō-Integral und das zweite Integral als ein gewöhnliches Riemann-Integral (über die stetigen Pfade des Integranden) zu verstehen.

Für den durch Y t = h ( W t ) {\displaystyle Y_{t}=h(W_{t})} für t ≥ ≥ 0 {\displaystyle t\geq 0} definierten Prozess lautet diese Darstellung in Differentialschreibweise

d Y t = h ′ ( W t ) d W t + 1 2 h ″ ( W t ) d t . {\displaystyle {\rm {d}}Y_{t}=h'(W_{t})\,{\rm {d}}W_{t}+{\frac {1}{2}}h''(W_{t})\,{\rm {d}}t\,.}

Ein stochastischer Prozess ( X t ) t ≥ ≥ 0 {\displaystyle (X_{t})_{t\geq 0}} heißt Itō-Prozess, falls

X t = X 0 + ∫ ∫ 0 t a s d s + ∫ ∫ 0 t b s d W s {\displaystyle X_{t}=X_{0}+\int _{0}^{t}a_{s}\,{\rm {d}}s+\int _{0}^{t}b_{s}\,{\rm {d}}W_{s}}

für zwei stochastische Prozesse a s {\displaystyle a_{s}} , b s {\displaystyle b_{s}} gilt (genaueres dazu unter stochastische Integration). In Differentialschreibweise:

d X t = a t d t + b t d W t . {\displaystyle {\rm {d}}X_{t}=a_{t}\,{\rm {d}}t+b_{t}\,{\rm {d}}W_{t}\,.}

Ist h : : R + × × R → → R {\displaystyle h\colon \mathbb {R} _{+}\times \mathbb {R} \to \mathbb {R} } eine in der ersten Komponente einmal und in der zweiten zweimal stetig differenzierbare Funktion, so ist auch der durch Y t := h ( t , X t ) {\displaystyle Y_{t}:=h(t,X_{t})} definierte Prozess ein Itō-Prozess, und es giltcite-ref-2[2]

d Y t = ∂ ∂ h ∂ ∂ t ( t , X t ) d t + ∂ ∂ h ∂ ∂ x ( t , X t ) d X t + 1 2 ∂ ∂ 2 h ∂ ∂ x 2 ( t , X t ) ( d X t ) 2 = ( ∂ ∂ h ∂ ∂ x ( t , X t ) a t + ∂ ∂ h ∂ ∂ t ( t , X t ) + 1 2 ∂ ∂ 2 h ∂ ∂ x 2 ( t , X t ) b t 2 ) d t + ∂ ∂ h ∂ ∂ x ( t , X t ) b t d W t . {\displaystyle {\begin{aligned}{\rm {d}}Y_{t}&={\frac {\partial h}{\partial t}}(t,X_{t})\,{\rm {d}}t+{\frac {\partial h}{\partial x}}(t,X_{t})\,{\rm {d}}X_{t}+{\frac {1}{2}}{\frac {\partial ^{2}h}{\partial x^{2}}}(t,X_{t})({\rm {d}}X_{t})^{2}\\&=\left({\frac {\partial h}{\partial x}}(t,X_{t})\,a_{t}+{\frac {\partial h}{\partial t}}(t,X_{t})+{\frac {1}{2}}{\frac {\partial ^{2}h}{\partial x^{2}}}(t,X_{t})\,b_{t}^{2}\right){\rm {d}}t+{\frac {\partial h}{\partial x}}(t,X_{t})\,b_{t}\,{\rm {d}}W_{t}\,.\end{aligned}}}

Hierbei bezeichnen ∂ ∂ h ∂ ∂ t {\displaystyle {\tfrac {\partial h}{\partial t}}} und ∂ ∂ h ∂ ∂ x {\displaystyle {\tfrac {\partial h}{\partial x}}} die partiellen Ableitungen der Funktion h {\displaystyle h} nach der ersten bzw. zweiten Variablen. Die zweite Darstellung folgt aus der ersten durch Einsetzen von ( d X t ) 2 = b t 2 d t {\displaystyle ({\rm {d}}X_{t})^{2}=b_{t}^{2}\,{\rm {d}}t} und Zusammenfassen der d t {\displaystyle {\rm {d}}t} - und d W t {\displaystyle {\rm {d}}W_{t}} -Terme.

Die Formel lässt sich auf n {\displaystyle n} Itō-Prozesse X = ( X 1 , … … , X n ) {\displaystyle X=(X_{1},\dots ,X_{n})} verallgemeinern. Sei h : [ 0 , ∞ ∞ ) × × R n {\displaystyle h:[0,\infty )\times \mathbb {R} ^{n}} in C 1 {\displaystyle C^{1}} in der ersten und C 2 {\displaystyle C^{2}} in den restlichen Variablen. Definiere Y ( t ) := h ( t , X ( t ) ) {\displaystyle Y(t):=h(t,X(t))} dann gilt

d Y ( t ) = ∂ ∂ h ∂ ∂ t ( t , X ( t ) ) d t + ∑ ∑ i = 1 n ∂ ∂ h ∂ ∂ i ( t , X ( t ) ) d X i ( t ) + 1 2 ∑ ∑ i , j = 1 n ∂ ∂ 2 h ∂ ∂ i ∂ ∂ j ( t , X ( t ) ) d [ X i , X j ] ( t ) . {\displaystyle {\rm {d}}Y(t)={\frac {\partial h}{\partial t}}(t,X(t))\,{\rm {d}}t+\sum \limits _{i=1}^{n}{\frac {\partial h}{\partial i}}(t,X(t))\,{\rm {d}}X_{i}(t)+{\frac {1}{2}}\sum \limits _{i,j=1}^{n}{\frac {\partial ^{2}h}{\partial i\partial j}}(t,X(t)){\rm {d}}[X_{i},X_{j}](t).}

Sei ( X t ) t ≥ ≥ 0 = ( X t 1 , … … , X t d ) t ≥ ≥ 0 {\displaystyle (X_{t})_{t\geq 0}=(X_{t}^{1},\dotsc ,X_{t}^{d})_{t\geq 0}} ein R d {\displaystyle \mathbb {R} ^{d}} -wertiges Semimartingal und sei F ∈ ∈ C 2 ( R d , R ) {\displaystyle F\in C^{2}(\mathbb {R} ^{d},\mathbb {R} )} . Dann ist ( F ( X t ) ) t ≥ ≥ 0 {\displaystyle (F(X_{t}))_{t\geq 0}} wieder ein Semimartingal und es gilt

F ( X t ) − − F ( X 0 ) = ∑ ∑ j = 1 d ∫ ∫ 0 + t ∂ ∂ F ∂ ∂ x j ( X s − − ) d X s j + 1 2 ∑ ∑ j , k = 1 d ∫ ∫ 0 + t ∂ ∂ 2 F ∂ ∂ x j ∂ ∂ x k ( X s − − ) d [ X j , X k ] s c + ∑ ∑ 0 < s ≤ ≤ t ( F ( X s ) − − F ( X s − − ) − − ∑ ∑ j = 1 d ∂ ∂ F ∂ ∂ x j ( X s − − ) Δ Δ X s j ) . {\displaystyle {\begin{aligned}F(X_{t})-F(X_{0})=&\sum _{j=1}^{d}\int _{0+}^{t}{\frac {\partial F}{\partial x^{j}}}(X_{s-}){\rm {d}}X_{s}^{j}+{\frac {1}{2}}\sum _{j,k=1}^{d}\int _{0+}^{t}{\frac {\partial ^{2}F}{\partial x^{j}\partial x^{k}}}(X_{s-}){\rm {d}}[X^{j},X^{k}]_{s}^{c}\\&{}+\sum _{0<s\leq t}\left(F(X_{s})-F(X_{s-})-\sum _{j=1}^{d}{\frac {\partial F}{\partial x^{j}}}(X_{s-})\Delta X_{s}^{j}\right).\end{aligned}}}

Hierbei bedeutet:

• X s − − = lim u ↑ ↑ s X u {\displaystyle \textstyle X_{s-}=\lim _{u\uparrow s}X_{u}} der linksseitige Grenzwert,

• das Integrationsgebiet 1 [ 0 + , t ] {\displaystyle 1_{[0+,t]}} bedeutet 1 ( 0 , t ] {\displaystyle 1_{(0,t]}} . Ein Semimartingal kann bei X 0 {\displaystyle X_{0}} einen Sprung haben, das heißt X 0 ≠ ≠ X 0 + {\displaystyle X_{0}\neq X_{0+}} und somit wird sichergestellt, dass nur über ( 0 , t ] {\displaystyle (0,t]} integriert wird und der Anfangswert F ( X 0 ) {\displaystyle F(X_{0})} wird deshalb nicht über das Integral gedeckt.

• Δ Δ X s j = X s j − − X s − − j {\displaystyle \Delta X_{s}^{j}=X_{s}^{j}-X_{s-}^{j}} der zugehörige Sprungprozess.

• Mit [ X j , X k ] c {\displaystyle [X^{j},X^{k}]^{c}} wird die quadratische Kovariation der stetigen Anteile der Komponenten X j {\displaystyle X^{j}} und X k {\displaystyle X^{k}} bezeichnet.

Falls X {\displaystyle X} ein stetiges Semimartingal ist, verschwindet die letzte große Klammer nach dem Plus und es gilt [ X j , X k ] c = [ X j , X k ] {\displaystyle [X^{j},X^{k}]^{c}=[X^{j},X^{k}]} .

Schreibt man den Ausdruck [ X j , X k ] t c := [ X j , X k ] t − − ∑ ∑ s ≤ ≤ t Δ Δ X s j Δ Δ X s k {\displaystyle [X^{j},X^{k}]_{t}^{c}:=[X^{j},X^{k}]_{t}-\sum \limits _{s\leq t}\Delta X_{s}^{j}\Delta X_{s}^{k}} aus, so erhält man für eine Funktion f ∈ ∈ C 2 ( R d , R ) {\displaystyle f\in C^{2}(\mathbb {R} ^{d},\mathbb {R} )} die Form

f ( X t ) − − f ( X 0 ) = ∑ ∑ j = 1 d ∫ ∫ 0 + t ∂ ∂ f ∂ ∂ x j ( X s − − ) d X s j + 1 2 ∑ ∑ j , k = 1 d ∫ ∫ 0 + t ∂ ∂ 2 f ∂ ∂ x j ∂ ∂ x k ( X s − − ) d [ X j , X k ] s + ∑ ∑ 0 < s ≤ ≤ t ( Δ Δ f ( X s ) − − ∑ ∑ j = 1 d ∂ ∂ f ∂ ∂ x j ( X s − − ) Δ Δ X s j − − 1 2 ∑ ∑ k , j = 1 d ∂ ∂ 2 f ∂ ∂ x j ∂ ∂ x k ( X s − − ) Δ Δ X s j Δ Δ X s k ) . {\displaystyle {\begin{aligned}f(X_{t})-f(X_{0})=&\sum _{j=1}^{d}\int _{0+}^{t}{\frac {\partial f}{\partial x^{j}}}(X_{s-}){\rm {d}}X_{s}^{j}+{\frac {1}{2}}\sum _{j,k=1}^{d}\int _{0+}^{t}{\frac {\partial ^{2}f}{\partial x^{j}\partial x^{k}}}(X_{s-}){\rm {d}}[X^{j},X^{k}]_{s}\\&{}+\sum _{0<s\leq t}\left(\Delta f(X_{s})-\sum _{j=1}^{d}{\frac {\partial f}{\partial x^{j}}}(X_{s-})\Delta X_{s}^{j}-{\frac {1}{2}}\sum _{k,j=1}^{d}{\frac {\partial ^{2}f}{\partial x^{j}\partial x^{k}}}(X_{s-})\Delta X_{s}^{j}\Delta X_{s}^{k}\right).\end{aligned}}}

wobei Δ Δ f ( X s ) := f ( X s ) − − f ( X s − − ) {\displaystyle \Delta f(X_{s}):=f(X_{s})-f(X_{s-})} .

→

Hauptartikel

:

Stratonowitsch-Integral#Itō-Formeln

Sei X = ( X 1 , … … , X n ) {\displaystyle X=(X^{1},\dots ,X^{n})} ein R n {\displaystyle \mathbb {R} ^{n}} -Semimartingal und f ∈ ∈ C 2 ( R n , R ) {\displaystyle f\in C^{2}(\mathbb {R} ^{n},\mathbb {R} )} , dann ist f ( X ) {\displaystyle f(X)} ein Semimartingal und es giltcite-ref-3[3]

f ( X t ) − − f ( X 0 ) = ∑ ∑ i = 1 n ∫ ∫ 0 + t ∂ ∂ f ∂ ∂ x i ( X s − − ) ∘ ∘ d X s i + ∑ ∑ 0 < s ≤ ≤ t ( f ( X s ) − − f ( X s − − ) − − ∑ ∑ i = 1 n ∂ ∂ f ∂ ∂ x i ( X s − − ) Δ Δ X s i ) . {\displaystyle f(X_{t})-f(X_{0})=\sum \limits _{i=1}^{n}\int _{0+}^{t}{\frac {\partial f}{\partial x_{i}}}(X_{s-})\circ dX_{s}^{i}+\sum \limits _{0<s\leq t}\left(f(X_{s})-f(X_{s-})-\sum \limits _{i=1}^{n}{\frac {\partial f}{\partial x_{i}}}(X_{s-})\Delta X_{s}^{i}\right).}

Hans Föllmer erweiterte die Formel von Itō auf (deterministische) Funktionen mit beschränkter quadratischer Variation.cite-ref-4[4]

Sei f ∈ ∈ C 2 {\displaystyle f\in C^{2}} eine reell-wertige Funktion und x : [ 0 , ∞ ∞ [ → → R {\displaystyle x:{[0,\infty [}\to \mathbb {R} } eine Càdlàg-Funktion mit endlicher quadratischer Variation. Dann gilt

f ( x t ) = f ( x 0 ) + ∫ ∫ 0 t f ′ ( x s − − ) d x s + 1 2 ∫ ∫ ] 0 , t ] f ″ ( x s − − ) d [ x , x ] s + ∑ ∑ 0 ≤ ≤ s ≤ ≤ t ( f ( x s ) − − f ( x s − − ) − − f ′ ( x s − − ) Δ Δ x s − − 1 2 f ″ ( x s − − ) ( Δ Δ x s ) 2 ) ) . {\displaystyle {\begin{aligned}f(x_{t})&=f(x_{0})+\int _{0}^{t}f'(x_{s-})\mathrm {d} x_{s}+{\frac {1}{2}}\int _{]0,t]}f''(x_{s-})d[x,x]_{s}\\&+\sum _{0\leq s\leq t}\left(f(x_{s})-f(x_{s-})-f'(x_{s-})\Delta x_{s}-{\frac {1}{2}}f''(x_{s-})(\Delta x_{s})^{2})\right).\end{aligned}}}

• Für Y t = sin ⁡ ⁡ ( W t ) {\displaystyle Y_{t}=\sin(W_{t})} gilt d Y t = cos ⁡ ⁡ ( W t ) d W t − − 1 2 sin ⁡ ⁡ ( W t ) d t {\displaystyle {\rm {d}}Y_{t}=\cos(W_{t})\,{\rm {d}}W_{t}-{\tfrac {1}{2}}\sin(W_{t})\,{\rm {d}}t} .

• Mit Hilfe der Formel kann man einfach beweisen, dass die geometrische brownsche Bewegung

S t = S 0 e r t − − 1 2 σ σ 2 t + σ σ W t {\displaystyle S_{t}=S_{0}e^{rt-{\frac {1}{2}}\sigma ^{2}t+\sigma W_{t}}}

eine Lösung der stochastischen Differentialgleichung von Black und Scholes d S t = r S t d t + σ σ S t d W t {\displaystyle {\rm {d}}S_{t}=rS_{t}\,{\rm {d}}t+\sigma S_{t}\,{\rm {d}}W_{t}}

ist.

Hierzu wählt man X t = W t {\displaystyle X_{t}=W_{t}} , also a t = 0 , b t = 1 {\displaystyle a_{t}=0,\;b_{t}=1} .

Dann ergibt die Formel mit h ( t , x ) = S 0 e r t − − 1 2 σ σ 2 t + σ σ x {\displaystyle h(t,x)=S_{0}e^{rt-{\frac {1}{2}}\sigma ^{2}t+\sigma x}} : d S t = [ ( r − − σ σ 2 2 + σ σ 2 2 ) S 0 e r t − − 1 2 σ σ 2 t + σ σ W t ] d t + [ σ σ S 0 e r t − − 1 2 σ σ 2 t + σ σ W t ] d W t = r S t d t + σ σ S t d W t . {\displaystyle {\rm {d}}S_{t}=\left[\left(r-{\frac {\sigma ^{2}}{2}}+{\frac {\sigma ^{2}}{2}}\right)S_{0}e^{rt-{\frac {1}{2}}\sigma ^{2}t+\sigma W_{t}}\right]{\rm {d}}t+\left[\sigma S_{0}e^{rt-{\frac {1}{2}}\sigma ^{2}t+\sigma W_{t}}\right]{\rm {d}}W_{t}=rS_{t}\,{\rm {d}}t+\sigma S_{t}\,{\rm {d}}W_{t}\,.}

• Ist ( W t ) t ≥ ≥ 0 {\displaystyle (\mathbf {W} _{t})_{t\geq 0}} ein d {\displaystyle d} -dimensionaler Wiener-Prozess und F : : R d → → R {\displaystyle F\colon \mathbb {R} ^{d}\to \mathbb {R} } zweimal stetig differenzierbar, dann gilt für Y t = F ( W t ) {\displaystyle Y_{t}=F(\mathbf {W} _{t})}

d Y t = ∇ ∇ F ( W t ) T ⋅ ⋅ d W t + 1 2 Δ Δ F ( W t ) d t {\displaystyle \mathrm {d} Y_{t}=\nabla F(\mathbf {W} _{t})^{\mathsf {T}}\cdot \mathrm {d} \mathbf {W} _{t}+{\frac {1}{2}}\Delta F(\mathbf {W} _{t})\,\mathrm {d} t} ,

wobei ∇ ∇ F {\displaystyle \nabla F} den Gradienten und Δ Δ F {\displaystyle \Delta F} den Laplace-Operator von F {\displaystyle F} bezeichnen.

Es gibt verschiedene Varianten von Itō-Formeln für unendlich-dimensionale Räume (z. B. Pardouxcite-ref-5[5], Gyöngy-Krylowcite-ref-6[6], Brzezniak-van Neerven-Veraar-Weiscite-ref-7[7]).

• Euler-Maruyama-Verfahren

• Philip E. Protter: Stochastic Integration and Differential Equations (2nd edition), Springer, 2004, ISBN 3-540-00313-4.

cite-note-11. ↑ Kiyoshi Itô: On a formula concerning stochastic differentials. In: Nagoya Math. J. Band 3, 1951, S. 55–65 (projecteuclid.org).

cite-note-22. ↑ Hui-Hsiung Kuo: Introduction to Stochastic Integration. Springer, 2006, ISBN 978-0387-28720-1, S. 103 (eingeschränkte Vorschau in der Google-Buchsuche).

cite-note-33. ↑ Philip E. Protter: Stochastic Integration and Differential Equations. Hrsg.: Springer. 2004, ISBN 3-540-00313-4, S. 277–278.

cite-note-44. ↑ Hans Föllmer: Calcul d'Ito sans probabilités. In: Séminaire de probabilités de Strasbourg. Band 15, 1981, S. 143–144 (numdam.org).

cite-note-55. ↑ E. Pardoux, E: Équations aux dérivées partielles stochastiques de type monotone. In: Séminaire Jean Leray. Nr. 3, 1974 (numdam.org).

cite-note-66. ↑ I. Gyöngy und N. V. Krylov: Ito formula in banach spaces. In: Springer, Berlin, Heidelberg (Hrsg.): Arató, M., Vermes, D., Balakrishnan, A.V. (eds) Stochastic Differential Systems. Band 36, 1981, doi:10.1007/BFb0006409.

cite-note-77. ↑ Z. Brzezniak, J. M. A. M. van Neerven, M. C. Veraar und L. Weis: Ito's formula in UMD Banach spaces and regularity of solutions of the Zakai equation. 2008.